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Abstract. We provide additional methods for the evaluation of the integral 

„ , f°° dx 

JVn 4 a;m := / r^- 

' V Jo (x 4 + 2ax 2 + l) m+1 

where m £ N and a £ (—1, oo) in the form 

7T 

No4(a:m) = -7- TTrrPmia) 

' K ' ' 2 m + 3 / 2 (a + l) m + 1 /2 mK ' 

where P m (a) is a polynomial in a. The first one is based on a method of 
Schwinger to evaluate integrals appearing in Feynman diagrams, the second 
one is a byproduct of an expression for a rational integral in terms of Schur 
functions. Finally, the third proof, is obtained from an integral representation 
involving modified Bcssel functions. 



1. Introduction 

The definite integral 

dx 



(1.1) N 0A {a-m) 
with m e N and a > — 1 has the value 



{x A + 2ax 2 + l) m +! ' 



( L2 ) N 0A(a; m) = 2m+3/2(g ^ 1)m+1/2 i 3 m(o) 

where 

(1 .3, P^ 2 -g 2 f :: f)(^) (o+1 > 



k 



The reader will find in [3] a survey of the many proofs of (|1.2[) that have appeared 
in the literature. After this survey was written, other proofs have appeared; see 
[51 IH1 [HI OH]- The goal of this note is to present some additional proofs: the first 
one is based in the Schwinger parametrization, a second one using Schur functions 
and the last one involves a representation involving integrals of Bessel functions. 
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2. SCHWINGER PARAMETRIZATION 

The method of Schwinger parameters, described in classical quantum theory 
books such as [TS] , is utilized here to give a proof of (|1.2p . A preliminary discussion 
is presented first. 

Starting from 

AiA 2 J Q [ X A 1 + [l - x) A 2 f Jo Jo (x 1 A 1 +x 2 A 2 ) 2 
where <5 is the Dirac measure, it can be deduced by induction that 

( 2 - 2 ) rr L ^ = {n -^[ ife^j^ n* 8 '- 

The induction step uses the n-times differentiated version of (|2.1[) , namely 

(2 - 3) a^ = LL { x lAl+ x 2 A 2 r^ dxldX2 - 

Repeated differentiations of (|2.2p yields the more general result 

i v{v) r s(i - Er=i^)nr=i^ -1 ^ 



(2 ' 4) nti at nr=i r m y [0 ,i]» (e? =1 m< 

with z/ = v\ + • • • + v n and A; G C. 

An alternative proof has been provided by Shapiro in his lecture notes [17] . This 
direct proof of (|2.2|) starts with 



(2.5) i = / e~ Ax dx 



o 



for RcA > 0. It follows that 



(2.6) IlT= < >: ' '" H^- 

Now let x = x\ H x n and j/j = to obtain 

n n / n \ 

(2.7) JJ dx, = x^dx Y[s(l-Y,y l ) d Vi 

i=l i=l \ i=l / 

so that 

n r n / n \ 

(2-8) n x = / n m 1 - e ^ N" -1 ^ e_x ELi Aiyi 

i=1 A ->R? i= i V i=l / 
After a change of variable, the latter integral is evaluated as 

,oo / n \" n 

(2.9) ^ a^e"* E "=* = r(n) f ^ A l2/l J 

resulting in 

i=1 ^ ^ i=1 12^=1 A^j 

Shapiro [17] states: 
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An interesting anecdote of physics history is that Schwinger re- 
mained bitter that a virtually identical mathematical trick be- 
came commonly known as Feynman parameters. Why two brilliant 
physicists, each of whom had an appropriately won a Nobel prize, 
should fight over what is essentially a trivial mathematical trick, 
is an interesting question in the sociology of physicists. 

An additional remark is that this representation is also known in probability the- 
ory, and was apparently introduced by Mauldon [13] . The n-dimensional Dirichlet 
distribution with parameters {vi} l<i<n of a random vector X reads 



f x (xi,...,x n )= Yi^r\vi) II ^' 



Xi 

i=l 



with v = V\ + . . . + v n . Mauldon proved that a Dirichlet distributed random vector 
X satisfies the following equality 

/ n \~ u n 

\i=l / i=l 

which is exactly (I2.4j) . 

This procedure is now applied to the integral No^(a;m) written as 

f°° dx 

(2.12) No4(a;m)= / — ^— ^- — — — - — — - dx 

y ' ' y ' J (1 + e i6 x 2 ) m+1 (1 + e- l9 x 2 ) m+1 

with a — cosO and 9 G (— ir, n). The corresponding parameters are n = 2 and 
V\ = v 2 = m + 1 . 

Recall the expression for the beta function 

(2.13) B(a,b)= [ x a - l {l-x) b - 1 dx 

Jo 

and its relation to the gamma function 

(2.14) B(a,b) T{a + b) 



T(a)T(b) 
to write 

»r r 1 f°° , f 1 x^{l- Xl ) m d Xl 

No a a; to = — / dx I ^ — nr- 

, V B(m + l,m + l) J J [ 1 + X 2 x ( Xl e*e + (i - Xl ) e -^)] 2m+2 

Formula 3.194.3 in [9] gives 

valid for u 6 C with |Arg(u)| < 7r. This yields 

B(|,2m+§) f 1 xf(l-xi) m dxi 



(2.16) N Aa\m)= , 

2B(m + l,TO + l) J ^f Xie i& + (i _ Xl ) e -iB 

The next lemma provides an evaluation of the integral in <|2.16j> . 
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Lemma 2.1. Define 

1 x m (l - x) m dx 



(2.17) I m {z) 
Then 



o y/l + (z 2 -l)5 



2k jm—h 



2 2m + 2 ^ (2m-2k\ fm + k\ 2k 

m(z)_ (-+i)r m +i )CD(-+i) 2m+i froV--fc ){ k r +l) z 

Proof. Expand the integrand in powers of x and integrate. The result is simplified 
using the value for the beta function 

a + b / a + b\ 1 
(2.18) B (a,b) = — 

for a, b eN. It follows that 



^(z) = V , 2 Hz 2 ~l) k x m+k {l-x) m dx 



( fc 2^B(m + fc + l,m + l)(z 2 -1) A 

l — n V * 



k=0 



— Y 

+ i f-r 



OO /r\;\ /n i 7 i i\ —1 /i 2 \ & 



2fc\ /2m + fc + l\ fl-z 



, fc / V m + 1 

fc=0 x 



The result now follows from the identity 

'2fc\ /2m + k + 1\ _1 fl-z' 



\ m + 1 



fe=0 

2 2m + 2 A (2m - 2k\ fm + k\ 2k m _ k 

CZt 1 ) CD (* + D 2m+1 fro V m - k) { k ) * + } z ' 

To establish (|2.19[) . the Wilf-Zeilberger automated method [16] is employed. It 
provides the recurrence 

(z 2 - l) 2 (4m + 9) (4m + 7)a m+2 - (z 4 - 6z 2 + l)(2m + 6) (2m + 3)a m+ i~ 

4z 2 (m + 2)(m + 3)a m = 

that both sides of (|2.19p satisfy. Checking the initial values at m = and m = 1 is 
elementary and the proof is complete. □ 

The evaluation of (|1.2p is restated in a slightly different form, as the content of 
the next statement. 



Theorem 2.2. Let m <G N, z = e ie and a = cos#. Then 

B{\,2m+l) 
2B{m + l,m + 1 



(2-20) N 0A (a; m) = + 2) 



holds. 
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3. A CONNECTION WITH SCHUR FUNCTIONS 

The integral 

(3.1) Gn(q) = -/ n^CT^ 

" , ° fc=l + q k 

is clearly a symmetric function of the parameters qi,...,q n . It has been expressed 
in [2] as 

(3-2) G„(q) = 7"- 1)( M ' 

e«(q) s A („)(q) 

where A(n) is the partition A(n) = (n — 1, n — 2, . . . , 1), e„(q) = gi^2 • • • Qn and the 
Schur function corresponding to a partition fi defined by 

a M +A(n)(q) 



«A(n)(q) 



(3.3) a M (q) 

with <z M (q) := det(gf j )i<ij< n - 

In the current problem, let a = \{w 2 + w~ 2 ) so that 

f°° dx 

(3.4) iVo4(o;m)= / -r-r. — , , , „ ^ — rr- 

Now take the partition A as above and the (2m + 2)— tuple q given by 

(3.5) A(n) = (n — 1, n — 2, . . . , 1), q = (w, . . . , w, u> _1 , . . . , w^ 1 ) 

where q has m + 1 copies of it) and u> . The expression p.2[) . that appears as 
Theorem 5.1 of [2], provides the next result. 

Theorem 3.1. The quartic integral is given by 

/, m , r i \ 7T s A(2m+l)(q) 

(3.6) N 0A (a;m) = - — . 

2 SA(2m+2)(q) 

As a corollary, this implies a specialization for the Schur functions. 
Corollary 3.2. If q = (w, . . . , w, w" 1 , . . . , w^ 1 ) is of length 2m + 2, then 
s\(2m+i){w, ...,w, w' 1 ,. . . , w- 1 ) _ 2 _ 2m /2fc\ /2m - fc\ 1 

Theorem 13.11 has a natural generalization to integrals of the form 

dx 



Q2n{x) m 1 

where 

(3.7) Q 2n {x) = x 2n + a lX 2n - 2 + a 2 x 2n - 4 + ■■■ + a 2 x 4 + a Y x 2 + 1 
is a palindromic polynomial. These polynomials factor as 

(3.8) Q2n(x) = {X 2 + wl){x 2 + Wf 2 ) • • • {x 2 + U> 2 )(x 2 + W~ 2 ). 

Consider the partition A and the n-concatenation q = (qi, ■ ■ ■ ,q n )j of the 2m- 
tuplcs qj, by 

(3.9) A(s) = (s-l,s-2,...,l), and qj = (wj, . . . , wf,wj l , . . . , wj x ) 
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where qj has m copies of each Wj and wj . A direct application of Q3.2p gives the 
next result. 

Proposition 3.3. Preserving the notation from above, 

(3 1Q^ S dx - 77 S A(2rcm-l)(<7) 



Q2n{x) m 2 Sx(2nm){q) 
4. A CONNECTION WITH BESSEL FUNCTIONS 

The quartic integral (ll.2[) is evaluated in this section using the modified Bessel 
function K a (x). The classical Bessel function J a (x) is usually defined via the 
differential equation 

(4 1) d 2 y^ldy f +°^\ 

dx 2 x dx \ x 1 ) 

The modified function I a (x) = e^ a7rl ^ 2 J a (xe 7Tl ^ 2 ) is also real valued. In the case 
a ^ Z, the functions I a (x) and I- a (x) form a basis for the solutions of (|4.ip . To 
deal with the situation of integer parameters, it is convenient to define 

(4.2) K a (x) = — [I_ a (x) - I a (x)] . 

I Sin7TQ 

This is the modified Bessel function of second kind, also called MacDonald function. 

The proof of (|1.2|) proceeds along the following lines: the first step is to obtain 
an expression for iVo,4(a, m) as an integral involving Kiu(x). An entry in |12j gives 
a hypergeometric form of this integral. The final step is to prove a hypergeometric 
identity that transforms this form to a result in [7] . Details about Bessel functions 
can be found in [4] and chapter 10 of |14) . 

The identity 

(4.3) fe-V^^, 

JO c 

valid for Re c > 0, is used with c = x A + 2ax 2 + 1 to produce 

(A 4) = / e -u u m -u(x A +2ax 2 ) j 

[ ' (x i + 2ax 2 + l) m + 1 r(m+l)7 

Integration yields 

-| pOO poo 

(4.5) N Qi {a-m) = — -/ e~ u u m e~ u ^ i+2ax ^ dudx. 

r(m + i)y J 

The restriction Re c > is satisfied by taking — 1 < a < 1. 
Entry 3.469.1 in |9] gives 

The choice [i — u and v = au yields: 

Lemma 4.1. The quartic integral No^aim) is given by 

l 

cym— ^ t>oo 

' —bt \ 



(4.7) A , 4 (a;m) = 3-/ t m e' bt K 1/4 (t) dt 



T(m + l)a 2m+ 2 



where b — 2/a 2 — 1. 
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This expression can be evaluated using |T21 vol.2, 2.16.6.2] 



(p + c) 



a + v, v + 



p — c 



p + c 



valid for Re (c + p) > 0, Re a > | Re^|. This yields the expression 



(4.8) N . A {a;m) 



/^a r(m+|)r(m+|) 



r (m - 



1*1 



f m + 


5 3 




4,4 




\ 771 


4- - 

T 2 





Using [1, 15.3.24] this can be written as 



(4.9) 



Nqa (a; m) 



0Fr(2777 +|) 

~2~ r (2m + 2) 



2*1 



2m 



3 1 

2' 2 
3 



1 - a 



-m, m + 1 

™ 1 3 



1 -a 



The proof of (|1.2[) is now based on the representation 

1 

(4.10) N QA (a ]m )= (a + i)m+1/2 2?(2m+|,l) 2 *i 

described in [7. In particular, the polynomial P m (a) is identified there as the Jacobi 
polynomial .P™ ' \ a )- 

Matching both representations for NQ A {a\m) shows that (|1.2I) follows from the 
next result. 



Proposition 4.2. The identity 



-m, m + 1 

™ _L 3 



(l-2z)(l-z) m+ 2 2 F X 



5 3 

, 4 3 4 



(4.11) 2F1 
holds. 

Proof. The WZ-method shows that both sides satisfy the recurrence 

(4.12) (4m + 7)(4m + 9)za m+2 + (2m + 3)(2m + 5)(4z 2 - 4z - l)a m+ i 

- (2m + 3)(2777 + 5)(z- l)a I( 



4z(l - z) 



0. 



To complete the proof, it suffices to check two initial values. The required iden- 
tities for m = and m = 1 are written in terms of t = 4z(l — z) in the form 



(4.13) 

and 

(4.14) 



2-^1 



5 3 
4' 4 

3 

2 



t = 



V2 



2F1 



9 3 

4' 4 
5 



(l-i) 1 /* Vl + Vl-t 
2^(3 + 2 V / T^t) 

5(1 + X /T _ 1)3/2^1— 

lementary using the fact 

t(l - %" - [c - (a + b + l)t] y - aby = 0. 
lalytic at t = with initial value 

t(l - + (| - 3t) y' - i|y = 



The proofs of these identities are elementary using the fact that the hypergeo- 
metric differential equation 

(4.15) 

has a unique solution analytic at t = with initial value y(0) — 1. Indeed, both 
sides of (|4. 13[) satisfy 

(4.16) 
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with value y(0) = 1. This establishes (j47T3|) . The same holds for (|4.14|) using 



(4.17) t(l - t)y" + (| - At) y' - %y = 



□ 

Note 4.3. The identity in Proposition l4.2l can also be established using C. Koutschan 
package HolonomicFunctions [TT]. Denote the left hand side of (|4. 1 1|) by Fl and 
the right hand side by F2. The command 

Annihilator[Fl, 5 [to]] 

finds the operator 

(7+4TO)(9+4TO)zS'^ + (3 + 2TO)(5 + 2TO)(-l-4z+4z 2 )5' m -(3+2TO)(5+2TO)(-l + z) 

with S m being the shift in the discrete parameter to; that is, S m g(m;x) = g(m + 
l;x). The package claims that this operator annihilates the left hand side of (|4.11l) . 
This is precisely the recurrence (|4.12l) obtained before. The command 

Annihilator[Fl, 5 [to]] == Annihilator[F2, 5 [to]] 
returns True, showing that the right hand side F2 satisfies the same recurrence. 
The initial conditions can also be determined automatically. 

5. A SECOND PROOF OF THE MAIN IDENTITY VIA BESSEL FUNCTIONS 

The previous section contains the identity (|4.ip giving the quartic integral No^(a, m) 
as an integral involving Bessel functions: 

l 

1 -bt , 



(5.1) iVo4(a;TO) = 3-/ t m e~ bt K 1/4 (t)dt. 

r(m + 1) a 2m+ 2 Jo 

The table of integrals [9] contains, as Entry 6.611.3, the special case to = 0: 

poo 

(5.2) J e- bt K lfi (t) dt = - 7f -^=— [(b + Vb^l) 1 ^ -(b+ Vb^l)- 1 ^ 
that can be written as 

(5.3) / e- bt K 1/4 (t)dt=- 



2 VT+a 1 

using b = a 2 /2 — 1. Differentiating to times with respect to 6, it follows that 

(-l) m 2 m - 1/2 d m fir a 3 / 2 
(5 - 4) ^ 4(Q ' m) = ml db^ I 2 

The proof of (II. 2|) is thus reduced to finding an analytic expression for the deriva- 
tives in (EI 



Proposition 5.1. There exists a polynomial Q m (a) such that 
d m ( a 3 / 2 \ _ ( — 1 ) m to ! a 2 "'+3/ 2 



<% m V\/l+~a/ 22m (l + a) m+1 /2 

Proof. The case to = holds with Q (m) = 1. To complete the inductive step, it 
is shown that if Q m {a) is a polynomial then the function Q m +\{a) defined by the 
relation 
(5.6) 

d ( {-l) m m\ a 2m +3/2 ^ \ _ (_l)m+l( TO+ 1)| a 2m+7/2 

^ ^ 22™ (l + a)™+i/2 Qm(a) J - (l + a)™+ 3 / 2 ° m+l(a) ' 
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is also a polynomial. 

The chain rule and ^§ = — \o? show that (15.61) is equivalent to 

(5.7) 2(m + l)Q m+ i(a) - [2(m + l)(a + 1) + (2m + 1)] Q ro (a)+ 

[2(a+l) 2 -2( a + l)] Q^(o). 

Thus Q m +i(a) is a polynomial in a. □ 



The proof of (|1.2[) is now reduced to checking that the polynomial Q m in the 
previous lemma is given by P m (a) defined in p. 31) . 

Theorem 5.2. The polynomial Q m (a) is the same as P m {a). 

Proof. It suffices to check that P m (a) satisfies the same recurrence as Q m (a). To 
this end, compare the coefficients of (1 + a) k on both sides of (I5.7[) . The result is 
equivalent to 

nk . iN /2m-2fc + 2\ /m + fc + l\ , , (2m - 2k + 2\ (m + k - 1 

2 fe (m + l) = 2 /i+1 (m + l) 



m-fc + 1 y \ m + 1 / \ m — k + 1 

'2m — 2k\ fm + k 
m — k ) \ m 

>fc+i / " > " 



2 fc+1 (2m + l) 
2 fe+1 (fc- 1) 
2 k+2 k 



2m-2k + 2\(m + k~l 
m — k + 1 J V m 



2m — 2fc\ / m + k 
m — k J \ m 



This however can routinely be verified. Simply divide through by ( 2 ™_^' c ) i™^**) 
and the statement reduces to a simple polynomial identity. □ 



A small variation of this proof of (|1.2|) is obtained by differentiating 

poo w 

(5.8) / e~ bt Ki(t)dt 



4 



2V4(6 + 1)V2 v /72+7F+i 



directly with respect to the parameter b. The first few examples suggest the next 
result. 

Lemma 5.3. There are polynomials S m , T m such that 
(5.9) ' 



db m V (6 +l)i/ 2 y/y/2 + y/b+l. 



, )m S m (b) + Vb+lT m (b) 

[ ' 22™(6+l) m +i/ 2 ( v / 2 + \/feTT) m +i/ 2 ' 

Proof. The proof is by induction on m, and is obtained upon differentiating the 
stated expression for the m-th derivative. The base case m = is obvious. Assume 
(|5.9|) holds for m. Then differentiating the right hand side of (|5.9|) generates the 
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recurrence 

(5.10) S m+1 (b) = 2V2(2m + l)S m (b) 



-(b + 1) [iV2S' m (b) - (1 + 6m)T m (6) + 4(6 + 1)1^(6) 
T m +i(b) = 3(2m+l)5 , m (6)+4 V / 2mT m (&) 
-4(6+l)(s^(6) + v^C(6)), 

where 6 = 2/a 2 — 1. The initial condition So (6) — 1 an d = yield the 



result. 

Lemma 5.4. Let m' = [m/2\. Define 



□ 



U m {b) 

and 
V m (b) 



(i + &)■■ 



- W 

V i. t 

k=0 



2m - 2k\ fm + *\ 1< *^ /21 f k 



g U'ii 12,(1 +4r 



Then 

(5.11) P m (o) = £/„,(&) + oV m (6). 

Proof. The polynomial P m (a) is decomposed into its even and odd part using 

(5.12) (1 + af = \ [(1 + af + (1 - a) k ] + i [(1 + a) fe - (1 - a) fc ] . 
The result follows. 



□ 



The proof of (|1.2j) now reduces to expressing the polynomials S m (&) and T m (b) 
in terms of U m (b) and This is given by 



(5.13) 
and 
(5.14) 



T m (b) = ^-J Um{b) 



if m is odd 



\/2a m ^ 1 I aV^(6) if to is even 



_ „, to!2 3 ™/ 2 \aV m (b) if to is odd 
S m (b) = — x > 



\f2a m [ U m (b) if to is even. 

The details are elementary and are left to the reader. 
Note 5.5. The identity (|jT3]l yields 



(5.15) 



N A{a,m) = 



TO !25m/2+3/2 ( fl _)_ ]\ro+l/2 



V2T m {b) + S m {b) 



THE EVALUATION OF A QUARTIC INTEGRAL VIA SCHWINGER, SCHUR AND BESSEL11 



6. Conclusions 
The value of the definite integral 

N ^;m) = J o {xi + 2ax2 + 1)m+ i> 

is given by 



where 




Several proofs of this identity appeared in the literature. The current work 
includes proofs relating this identity to Bessel functions, Schur polynomials and a 
method of Schwinger for the evaluation of definite integrals. 
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